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Abstract 

We study the effective interaction between differently charged polyelectrolyte-cohoid complexes 
in electrolyte solutions via Monte Carlo simulations. These complexes are formed when short and 
flexible polyelectrolyte chains adsorb onto oppositely charged colloidal spheres, dispersed in an 
electrolyte solution. In our simulations the bending energy between adjacent monomers is small 
compared to the electrostatic energy, and the chains, once adsorbed, do not exchange with the 
solution, although they rearrange on the particles surface to accomodate further adsorbing chains 
or due to the electrostatic interaction with neighbor complexes. Rather unexpectedly, when two 
interacting particles approach each others, the rearrangement of the surface charge distribution 
invariably produces anti-parallel dipolar doublets, that invert their orientation at the isoelectric 
point. These findings clearly rule out a contribution of dipole-dipole interactions to the observed 
attractive interaction between the complexes, pointing out that such suspensions can not be consid- 
ered dipolar fluids. On varying the ionic strength of the electrolyte, we find that a screening length 
K~^, short compared with the size of the colloidal particles, is required in order to observe the at- 
traction between like charged complexes due to the non-uniform distribution of the electric charge 
on their surface ('patch attraction'). On the other hand, by changing the polyelectrolyte/particle 
charge ratio, ^s, the interaction between like-charged polyelectrolyte-decorated (pd) particles, at 
short separations, evolves from purely repulsive to strongly attractive. Hence, the effective inter- 
action between the complexes is characterized by a potential barrier, whose height depends on the 
net charge and on the non-uniformity of their surface charge distribution. 



I. INTRODUCTION 



Understanding complexation of linear polyelectrolyte chains with oppositely charged 
particles has high potential for applications in nanoscience and bio-technology 2|-|7|. In a 
number of technological processes and biological systems, linear polyelectrolytes associate 
with a range of different nanoscopic particles or supramolecular assemblies, including 
inorganic particles, protein assemblies, micelles and liposomes, forming in aqueous media 
charged colloidal suspensions. DNA, which is a polyelectrolyte, provides a typical example 
of the relevance of these interactions. Within the cell nucleus the very long DNA chains 
are packed and ordered into much smaller structural units (nucleosomes) due to their 
interaction with strongly alkaline proteins, the histones, acting as 'spools' around which 
DNA winds. 

More in general, the long-range character of the electrostatic interaction, and its interplay, 
in a polyelectrolyte solution, with entropic and conformational effects, result in a very 
complex and interesting phenomenology 8|, |9|. Particularly, when polyelectrolytes adsorb 
on an oppositely charged extended surface, entropic, geometric and correlation effects play 
a pivotal role. 

There is now accumulating evidence that the polyelectrolyte adsorption^ due to the 

repulsion between the like-charged chains, occurs in a correlated manner lOJ, and that 



the intriguing phenomenon of the overcharging [Ill4l4|. when more polyelectrolyte than 
would be needed to neutralize the surface adsorbs, so that the sign of the net charge of the 
'decorated' surface is inverted, can be understood in terms of correlated adsorption [13| and 
'charge fractionalization' 15|, |l6|. The basic idea of 'fractionalization' is that by forming 
loops or leaving dangling ends, the adsorbed chains gain some conformational entropy. The 
vacancies left by these defects can locally be large enough to drive closer to the surface 
oncoming polyelectrolytes, that may adsorb, again with some loop and dangling ends (above 
the surface). As a result, instead of having a Z-charged chain in solution and the surface 
neutralized by a uniform coating of ordered chains laying fiat on it, Z 'disconnected' charges 
appear protruding from the adsorbed polyelectrolyte layer: the charge of the polyelectrolyte 
is fractionalized along the surface. This is an energetically favored configuration for 
the system, that gains some conformational entropy of the adsorbed chains and also, 
by 'diluting' its charge on the whole surface, part of the electrostatic self-energy of the 



Z-charged polyelectrolyte. 

By increasing the polyelectrolyte/particle charge ratio, the net charge of the decorated 

particles changes progressively from the value of the bare particle to the maximum 

'overcharging' inverting its sign at the isoelectric point. 

The correlated adsorption of the chains also results in a non-uniform distribution of the 

electrostatic charge on the surface of the decorated particles. As a consequence of this 

non-uniformity, an attractive interaction can be observed between like-charged particles 

This attractive component of the inter-particle potential, and the progressive variation 
of the net charge of the decorated particles, down to the neutralization point and up 
again to the maximum overcharging, that modulate the electrostatic repulsion, combine 
together yielding a peculiar phenomenon of 're-entrant condensation' of the polyelectrolyte- 
decorated particles. As the polyelectrolyte/particle ratio is increased, associated to the 
progressive reduction of the net charge of the primary polyelectrolyte-decorated particles 
(pd-particles), larger and larger clusters are observed. Close to the isoelectric condition 
the aggregates reach their maximum size, while beyond this point any further increase of 
the polyelectrolyte-particle charge ratio causes the formation of aggregates whose size is 
progressively reduced. 

Eventually, when the surface of the particles is completely saturated by the adsorbed 
polyelectrolyte, i.e. the overcharging has reached its maximum extent, the size of the 
particles in the suspension equals again the size of the original colloidal particles, plus 
a thin layer of adsorbed polymer [^, |5|, ll9|. From here on, by further increasing the 
polymer /particle ratio, the excess polyelectrolyte, that does not adsorb any more on the 
particles, remains 'freely' dissolved, contributing to the overall ionic strength of the solution. 



This intriguing phenomenology has been observed in a variety of polyelectrolyte-col 
systems dispersed in aqueous solutions such as polyelectrolyte- micelle complexes. 



latex particles 
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22|, dendrimers 
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23| . ferric oxide particles |2j|, phospholipid vesicles 
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(liposomes) 25l-l27l|. 'hybrid niosome' vesicles 28 1. 



fl 



Velegol and Thwar [l| have recently suggested an analytical model (based on the 



Derjaguin approximation and on an extension of the Hogg-Healy-Fuerstenau model 29l |) 



for the potential of mean force (VT-potential hereafter) between non-uniformly charged 
surfaces (planar and spherical), showing that the interaction between non-uniformly charged 
surfaces results in an inter-particle potential that, even in the case of like-charged particles, 
has an attractive component. 

Accurate simulations on a system composed by spherical particles interacting via VT- 
potential 30|| showed that the overall phenomenology is, at least qualitatively, well-described 
by this model. In particular, the balance of attractive and repulsive electrostatic interactions 
yields a potential barrier whose height increases with the size of the aggregates, justifying 
the stabilization of the aggregates. Aggregates stop growing after reaching a characteristic 
size, determined by the net charge and the degree of non-uniformity. Such behavior is 
consistent with the observed 'reentrant condensation' and with the thermally activated 
character of the aggregation 30, l3l[ |. 



Most of the previous works on polyelectrolyte-colloid complexes focused on the complex- 
ation behavior, i.e on the adsorption of PE chains on a single oppositely charged colloidal 



particle 
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38| or viceversa of one long PE on many particles 35|, l36| , being mainly con 



cerned with the overcharging phenomenon. A few authors discussed the interactions between 



two macroions in the presence of one 



chain 



32 



39| or several short chains 40 



ong (compared to the particle size) polyelectrolyte 



41| in different environments. Here we mainly focus 



instead on the effective interaction between decorated particles, considering the complexes 
formed by one spherical macroion and several PE-chains whose length is small compared 
with the particle size. We investigate the effective interaction between two of such deco- 
rated colloids in and electrolyte at varying the screening length, k~^, and the polyelectrolyte 



concentration. In our simulations we reproduce the conditions that are typical 



br the 



polyelectrolyte-liposome systems employed in several experimental works |19l. l25l. l28l. |42|-|45| 
and that show an interestingpotential for biotechnological applications as drug carriers for 
intra-cellular drug delivery |2[. To this aim we perform Monte Carlo (MC) simulations of 
PE-decorated particles, for different polyelectrolyte/particle charge ratios, in the limit of 



bending energy smal 
Dzubiella et al. 



compared to the electrostatic interaction (rather flexible chains), 
showed that, for rather small colloidal particles (with a size not 



exceeding ten times the radius of a polyelectrolyte monomer) and almost neutral complexes, 
the shape of the interaction potential as a function of the inter-particle distance is only 



minimally affected by the ionic strength of the solution. On the contrary, here we show 
that, for short chains and small monomers and for non-neutral complexes, this parame- 
ter plays, through the screening length k, a key role, modulating the balance of repulsive 
(due to the net charge) and attractive (patch) interactions. We find that the addition of 
the polyelectrolyte chains on the macroions' surfaces switches the short distance interaction 
from repulsive to attractive, yielding a potential barrier modulated by the net charge and 
by the non- uniformity of its distribution. The nature of the short range attractive inter- 
action between the decorated particles will be further clarified by calculating their dipole 
moments. Rather unexpectedly, the dipole moments of approaching particles remain anti- 
parallel, ruling out a contribution of dipolar interactions to the attractive component of the 
inter-particle forces, and showing that the condensation of the complexes is only controlled 
by the charge mismatch between the two opposing surfaces. 

II. THE MODEL 

We perform MC-simulations adopting the scheme proposed by Kremer et al. |46|, |47| for 
polyelectrolyte chains, and the pearl necklace model to generate off-lattice 3-dimensional 
polymer chains represented as a succession of A^ = 20 freely jointed Lennard- Jones spheres 
(beads). Each bead is a physical monomer with one positive elementary charge positioned at 
its center, and a diameter equal to one Bjerrum length, Ib = e^/eksT, the distance between 
two point charges e where their electrostatic interaction energy in a medium with dielectric 
permittivity e reduces to the thermal energy ksT (7.14 A in water at 25 °C). 
The fraction of ionized monomers / is set to 1 whereas the bond length is a fluctuating quan- 
tity whose oscillation is tuned by the binding potential. Assuming good solvent conditions, 
a shifted Lennard- Jones (LJ) potential is used to describe the purely repulsive excluded 
volume interaction between the A^ monomers: 
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r > 2lR„^ 

where Rm is the radius of the beads (the microscopic length scale), r their relative dis- 
tance, and e^j sets the energy scale. The connectivity of the bonded monomers is assured by 



a finite extension nonlinear elastic (FENE) potential 39|] acting between neighboring beads 
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where kp denotes the spring constant. Following Dzubiella et al. 39|] we have chosen 



kp = 7.0eLj, the maximal relative dis plac ement of two neighboring beads Vmax = ^Rm and 



the L-J parameter ej ] = 1.2kBT J39|, |48|, |49|]. These values of the parameters prevent the 
chains from crossing 48 1. 



Finally, electrostatic interactions between monomers have been taken into account via 
the Debye-Hiickel potential 



VDHir) 
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where z=l is the monomer valence, e the elementary charge and e is the dielectric per- 
mittivity of the medium, set to 7.08 ■ 10~^° j[^. n~^ is the Debye screening length, related to 
the electrolyte concentration as follows 



1000 ■ c^Na Y^ 
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where where Na is the Avogadro's number, Zi the valence of the ionic species, Cj the elec- 
trolyte concentration in mol/1, and ksT the thermal energy. The solvent is only taken into 
account via the dielectric constant. It must be noted that in the absence of any added salt, 
the system contains the small counterions resulting from the ionization of both the chains 
and the colloidal particles which, as stated above, is assumed to be complete. However, also 
these small counterions are not considered explicitly in the simulations but only taken into 
account via the screening constant n~^. 

We do not add a bending term to the total pair interaction energy between adjacent 
monomers. 

Using these parameters, the bond-length fe, defined as the distance between the centers of 
two consecutive monomers, has a weak dependence on ionic strength {h varies from 0.84 nm 
to 0.82 nm on passing from k = 1 nm~^ to k = 0.1 nm~^ ). In this conditions, the value of 
the Manning condensation parameter Is/b 50|, l5l|] is less than 1, justifying the choice / = 1 
(complete ionization). Hence, the total valence of one PE is simply given by Zp = Nz. 



The interactions between each monomer and spherical macroions are modeled via Debye- 
Huckel and truncated LJ potentials as follows: 

_ Z2;e^exp[-/t(r-i?,)] 
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r>2-e{Rrn + Rc) 

where Re is the radius of spherical macroions. The first potential (Eq. |5]) drives the 

adsorption process and determines the internal structure of the PE-layer, the second one 

(Eq. |6]), where e^j is set to lOe^j, reproduces the strong excluded volume effects due to the 

steric interaction between the monomers and the macroion. 

In our simulations we used, for different purposes, spherical particles with radii Re = 3.41 

nm and Re = 5 nm, as will be shown in the next section. 

In an attempt to reproduce what occurs experimentally when some oppositely charged 



a 



polyelectrolyte is added to a colloidal aqueous suspension (2|], we started simulating the 
formation of a decorated particle. A colloidal particle was located at the center of a cubic 
box of side Lbox = 60 nm. Then Nc chains were randomly placed within the box, outside 
the excluded volume of the central sphere. The macroion was kept fixed, while the polymers 



where moved according to the standard MC Metropolis algorithm 52| . Simulations have been 
carried out for a long equilibration time, typically of the order of 10^ MC steps, where one MC 
step consists of a trial translational move of every monomer. Periodic boundary condition 
have been applied. In the equilibrium final configuration, Nc chains are adsorbed on the 
oppositely charged spheres. The charge ratio ^s = \fNNc/Z\, where Z the valence of the 
spherical macroion, provides a measure of the net charge of the complex. This preparation 
procedure has been carried out for each set of parameters used in the simulations. 

Once a single complex, a pd-particle, was formed, we duplicated the whole structure 
(bare colloid + adsorbed PE chains) simply by translation. Then, to calculate the interac- 



tion energy [39|, l53|], the system was equilibrated, keeping fixed the macroions, at various 
distances. The surface-to-surface distance between the original particle and the copy was 
varied from 50 nm to nm in the first series of simulations (for the longer screen distance, 
K'^= 10 nm), and from 10 nm to nm in a second series (k~^ < 2 nm). 

8 



The electrostatic interaction between the macroions was taken into account using the Debye- 



Hiickel potential [39 1 



while the excluded volume interaction was again modeled via a LJ potential, as in Eq. |6l 
where Rm has to be substituted by Re- All the simulations involving two pd-particles (sec- 
tions Hill and [IV]) were performed using a cubic box much larger than the typical 'interaction 
volume' defined by the couple of particles; thus we did not need to apply periodic boundary 
conditions. 

To favor the rearrangement of the chains toward the global energy minimum (at fixed 
macroion-macroion distance) we adopted two kinds of trial moves: 

• translational moves of single monomers, 

• orientational moves of the whole polyelectrolyte layer around the macroion center. 

The diffusion process along the axis connecting the macroions' centers depends on the net 
charge of each complex and on reciprocal distance, as we shall see studying the polarization 
effects. One orientational move is performed randomly and, on average, every 5 translational 
MC-steps. When such move is performed, each spherical PE-layer is rotated around the 
center of the respective adsorbing macroion. Both the translational and orientational moves 
were applied with an acceptance ratio of 50 %. 

The two-step procedure (formation of a decorated particle plus duplication) employed 
to obtain a couple of interacting complexes, can be justified considering that the diffusion 
of the colloidal particles and the adsorption process occur on two different time scales. 
The huge surface/volume ratio which characterize colloids, and the distribution of this 
surface within the whole volume of the host phase, dramatically speeds the adsorption 
process. Upon mixing the polyelectrolyte solution and the particle suspension, the time 
requested to the polyelectrolyte chains to reach by diffusion the particles' adsorbing surface 
is significantly smaller than the diffusion time of the macroions and hence the effective 
interaction is essentially between already fully decorated particles. This conjecture, which 
appears in itself reasonable, has been recently substantiated with experimental evidence 
by Volodkin et al.J5|. By using different mixing protocols, varying the agitation speed and 



the order of mixing, these authors showed that the polyelectrolyte adsorption is almost 
immediate when compared with the characteristic times of the aggregation process, that 
are typical of the particles' diffusion. 



III. POLYELECTROLYTE ADSORPTION: THE EFFECTS OF SALT AND 
CHARGE RATIO 

The adsorbed polyelectrolyte (PE) chains play in principle a double role in determining 
the overall mean-force potential between the decorated particles: by building up a non 
zero dipole moment they can contribute repulsive or attractive components, moreover, at 
a sufficient short distance they produce a soft repulsive interaction. Both these effects 
depend on the distribution and on the conformation of the adsorbed chains. Thus, to 
gain insight on how the PE-layer of a non-neutral complex reacts to the presence of an 
approaching pd-particle, for different screening lengths k~^, is particularly important. The 
Debye-Huckel potential, accounting for the screening created by the presence of small ions 
in the solution, greatly simplifies the study of the effect of changes of the ionic strength on 
the PE-layer thickness. 

The typical value for the macroion charge density in cationic liposomes is E/) = 1.67 



nm~^ 5J]. We have selected this value for our study, resulting in a macroion charge of 
Z = —244 for Re = 3.41 nm and Z = —523 for R^ = 5 nm. In particular, the effect on the 
PE-layer of the screening length is here investigated using macroions with Re = 3.41 nm, 
and a surface charge ratio equal to 0.57, equivalent to 7 chains on each macroion (see the 
snapshot in the upper panel of Figure [3]). 

We calculated the thickness (6) of the PE layer as the mean distance between the 
monomers and the macroion surface. For different values of k, a common behavior of 
the thickness (6) as function of the distance H between two complexes is recovered. 
Figure [1] (panel A) shows that the approach of two complexes induces a weak but definite 
deformation of the polyelectrolyte adsorbed layer. When the distance between the two 
complexes becomes smaller than ^ Anm, the adsorbed monomers begin to be attracted 
by the field generated by the approaching complex. The resulting increase of their mean 
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FIG. 1: PANEL A: Averaged mean distance between the monomers and the bare macroion surface 
as a function of the distance K between two approaching macroions. Simulations are carried out 
for PE-coUoid complexes with R^ = 3.41 nm, ^g = 0.57 and different values of k. Black filled points 
-^ K = 0.1 nm~^; circles -^ k = 0.5 nm~^; squares -^ k = 0.75 nm~^; triangles — )> k = 1 nm~^; 
lozenges -^ k = 1.25 nm~^. PANEL B: Averaged mean distance between the monomers and the 
bare macroion surface as a function of k for isolated complexes. Inset: Average mean distance (6) 
for nm~^< k < 0.3 nm~^ in a log-log scale. 

distance from the macroion surface, the layer 'swelling', depends on the inverse screening 
length K. But, as expected, the distance where this increase appears, does not depend on k 
(both the monomer-colloid adsorption energy and the attraction induced by the approaching 
complex scale with the same n). At much shorter distances {H < 1 nm) the layers on 
the approaching complexes become thinner, each layer being squeezed by its homologous 
on the other complex. The position of the maximum of {6{H)) can be considered the 
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effective range of the repulsive steric interaction between the complexes that, as will be 
shown in Fig. |3l is K-independent, at least in the range of Debye constant values investigated. 



For comparison, Fig. [T] (panel B) shows the variation of the adsorbed layer thickness (6) 
of isolated pd-particles. The dependence of (6) on n points out that for the chosen values 
of the surface charge density on the macroions, S, and of the valence of monomers z, the 
system is in the so called 'screening-reduced' regime, as predicted by previous analytical 



and numerical studies for electrostatically driven polymer adsorption |55l |56| . Our data are 
consistent with the scaling laws obtained by Borukhov et al. 57|] via mean field arguments 
for the dependence of the PE-layer thickness on salt concentration: (6) ~ k, for high values 
of K below the desorption threshold {ude — l-5nm~^), while for low salt concentration the 
thickness tends to be unperturbed by the increase of the screening, i.e., (6) ~ kP (see the 
inset of figure [1]- panel B). 



We also investigated the effect on the PE-layer thickness of the polyelectrolyte-particle 
charge ratio. In a series of simulations we fixed the radius of the macroions to Re = 5 nm 
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FIG. 2: Averaged monomer-colloid distance (panel A) and average end-to-end distance of chains 
(panel B) for different charge ratio S,s, k = 1 nm ~^, Re = 5 nm^^. Dashed lines are guide to the 
eye only. 

and the Debye constant to k = 1 nm~^, and added a different number of PE-chains so that 
the polyelectrolyte/macroion charge ratio varies from ^^ = 0.11 (3 chains) to C,s = 1-5 (39 
chains), where the polymer is in large excess. 
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The increased roughness of the PE-layer and the consequent increment of its thickness 
(see Fig. |2l panel A) for relatively high values of C,s is a direct consequence of the lateral 
repulsion between the chains. The strong non linear behavior of (6) is due to the presence of 
a disordered three-dimensional arrangement of the chains, that do not lay flat on the surface 
but adopt a conformation characterized by loops and dangling ends. 
To estimate a typical size of the adsorbed coils we calculate the averaged end-to-end distance 



R. 
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(8) 



where rj(0) and ri{N) are the positions of the first and last monomer along the ith chain, 
respectively. 

Figure [2] (panel B) shows that at very low polyelectrolyte content, due to the steric constraint 
contributed by the spherical macroion, and to the negligible reciprocal interaction, the chains 
assume a rather extended conformation. The addition of more chains, due to the increased 
inter-chain repulsions, results in a reduction of the size of the adsorbed coils. This effect is, 
of course, directly correlated to the stiffness of chains: the introduction of an elastic bending 
energy would reduce the effect of lateral interactions on the coil size, that becomes negligible 
in the limit of rigid polymers. At higher polyelectrolyte content the shrinking of the chains 
levels off and large loops and dangling tails protruding from the surface appear, after that 
the condition of excess polyelectrolyte is reached (chains do not adsorb any more). The 
passage from a 'flat adsorption' to a 3D one, results in an extension of the steric repulsion 
range, that determines the shift of the global minimum of the effective interaction towards 
higher values of the inter-particle distance H. 

IV. EFFECTIVE INTERACTION 

For each value of k, and $,s, we calculated the mean force potential between two complexes 
at a distance H by averaging the variation of the total energy E of the system when the 
distance between the particles is reduced from infinity to H 

^H) = {E{H) - Eioo)) (9) 

As a first step, we investigated the role of the screening on the interaction between non 
neutral complexes (the parameters are the same as in section UTTl with ^g = 0.57, Re = 3.41 
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FIG. 3: Mean force potential between two spherical PE-colloid complexes with Re = 3.41 nm 
(,s = 0.57 and k = 1 nm ~^. The dashed line represents the potential calculated from Eq. \W\ The 
snapshot shows two interacting pd-particles for ^s = 0.57, H = 5 nm and Re = 3.41 nm 



nm and 7 chains adsorbed on each macroion). 

For large enough screening lengths, i.e. when kR^ < 1, the decorated particles are able to 
'see' each other completely, in other words all the charges on the two complexes interact with 
each other, and all of them contribute to the effective pair potential. In these conditions the 
interaction between the pd-particles could be described as a "particle-particle" interaction. 
Figure [3] shows the behavior of the mean force potential as a function of the inter-particle 
distance in these conditions, i.e. for k, = 0.1 nm~^ {uRc = 0.34): the interaction is purely 
repulsive. The short range attractive component due to the matching between oppositely 
charge patches on the approaching particles does not balance the strong repulsion between 
the like-charged complexes. 

In fact, the spatial profile of the mean force potential obtained can be exactly superimposed 
on the Debye-Huckel potential 

{Z -^sZfe^exp[-KH] 



Vdh{H) 



(10) 



A-Ke{H + 2Re){l + KRef 
where Z — C,sZ is the net charge of the complexes (Figure [3]). As expected, for low screening 
conditions, i.e. k{Rc + 5) < 1, where 6 is the PE-layer thickness and 6/Rc ^ 1, the 
interaction is totally governed by the monopole-monopole repulsion. In this regime the 
non-uniformity of the charge distribution on the particles' surface has negligible effects and 
the interaction can be simply described as Debye-Huckel. 
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On the contrary, in high screening regime, due to the charge non- uniformity, a short 
range electrostatic attraction appears, which can not be justified within the simple Poisson 
Boltzmann theory for like charged colloids. 

We performed different MC simulations calculating the total mean- force potential ^{H) for 
increasing values of k. Figure H] shows the effective interaction ^{H)/kBT for 0.5 nm~^< 
K, < 1.25 nm~^. 
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FIG. 4: Total mean force potential between two spherical pd-particles with Re = 3.41 nm and 
^s = 0.57 for different values of k. (o) k = 0.5 nni~"^; (D) k = 0.75 nm~^; (A) k = 1 nni^^; (0) 
K = 1.25 nm~^. The dashed line corresponds to the Debye-Huckel potential (Eq. [T0|) shown in 
Fig. [3] for K = 0.1 nm~^. 



The observed behavior of the mean-force potential as a function of the ionic strength can 
be intuitively understood if the potential is explicitly thought of as the difference between 
the total energy of the two complexes when they are isolated E{oo) and when they are at a 
distance H (see Eq. |9]). For an unscreened potential {kRc <^ 1) each complex interacts with 
all the charged elements of the opposing complex, hence, the averaged interaction between 
the two particles results in a function of their overall net charge. On the contrary, when 
the screening is sufficiently large {uRc > 1), the effective interaction is in practice limited 
to those parts of the complexes that are sufficiently close each other [H < k~^ < Re). The 
interaction could be hence defined as a surface- surface interaction (as opposite to a particle- 
particle one). In these conditions, the local attraction between oppositely charged domains 
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on the two opposing particles can give rise to a net attraction of purely electrostatic nature, 
even though the particles are, as a whole, like charged. 

For the charge ratio employed here (^s = 0.57), the crossover from repulsive to attractive 
regime is observed for 2.56 < nRc < 3.41. Above this range a complete destabilization of 
the dispersion occurs. 

At very short distances the effective interaction 'feels' the effect of the steric repulsion 
due to the overlapping of opposing PE-layers. The soft shell of adsorbed polyelectroljd;es 
enveloping the macroion 'hard sphere', significantly contributes to the interaction only when 
the two complexes are in close contact. 



H 


200- 

0- 

-200- 

-400- 


q ^ 
^^"^ ^~ 

^ ^■AA-A'B-a=-.Q-.-; 

v2'^ I".:*''"" 
•, ' ^- ■ 

.'p'V 

' / ' 

,'-v , , , 


Bare colloids 




.---¥--- -V - * 








- 


9* ; 




H[nm] 





H[nm] 



FIG. 5: The electrostatic component of the nieanforce potential <^e between PE-coUoid complexes 
for K = 1 nm^^, Re = 5 nm, Z = —523 and different values of ^s- Main panel: Dashed line — >■ 
potential calculated from Eq. [TU] for bare Z-macroions {S,s = 0); o) S,s = 0.11; D) S,s = 0.26 ; A) 
^s = 0.38; V) ^s = 0.46; 0) 6 = 0.57; •) ^s = 1-0; ■) 6 = 1-2; T) Cs = 1-33. Inset: potential 
profiles in the 'intermediate' i/-region, where an energy barrier between the non neutral complexes 
arises. 

In practice, the effect of the adsorption of the PE-chains is double: the progressive 
neutralization of the net charge of the complexes that progressively quenches the long range 
repulsion, and (due to the correlated character of the adsorption and to the presence of 
a screening) the appearance of a purely electrostatic short range attraction. In Figure 
|5] we note the progressive reduction of the repulsion in the electrostatic component of the 
effective interaction $e(-f^)/^_B^ as the polyelectrolyte content is increased and, as expected. 
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the rising of a strong short range attraction for ^s > 0.46, well below the isoelectric point 

In the inset of Figure [5] three different profiles of the electrostatic interaction $e(-^)/^_B^ 
are shown, to emphasize what occurs below and above the isoelectric point {^^ = 1.0), where 
there is not monopole-monopole repulsion at all. When a patch attraction is present but 
the complexes are not neutral, the competition between the residual repulsive interaction 
and the shorter ranged attraction determines a potential barrier, whose strength depends 
on the balancing between these two components. The presence of such barrier justifies [30(] 



the aggregation behavior observed in several co 
amount of added polyelectrolyte chains [2, | l9 



loid-polyelectrolyte systems regulated by the 
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FIG. 6: Effective interaction between PE-colloid complexes for k = 1 nin ^, Re = 5 nm, Z = —523 
and different values of ^s- Main panel: o) ^^ = 0.11; D) ^^ = 0.26 ; A) ^^ = 0.38; V) ^s = 0.46; 0) 
^^ = 0.57; •) ^^ = 1.0; ■) ^, = 1.2; t) ^s = 1-33. Inset: potential profiles (6 = 0.57, 1.0, 1.2) in 
the region of the global minimum. 

By adding a 'steric' contribution to the inter-particle potential, due to the overlapping of 
PE-layers, we recover the complete effective interaction between the decorated macroions. 
Figure [6] shows different profiles of ^{H)/kBT for 0.11 < ^s ^ 1.33. The global minimum 
of the potential is the result of the competition between the repulsive steric component and 
the patch attraction. The position of this minimum is affected by the charge ratio only for 
high values of C,s, where the PE-layer cannot be considered fiat anymore and the dangling 
ends protruding from macroion increase the overlap volume of the adsorbed layers on the 
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opposing particles, thus causing a shift of the minimum towards higher values of H. 

V. THE DIPOLE MOMENT OF THE COMPLEXES: THE ANTI-PARALLEL 
DOUBLETS 

As a further investigation of the effects of the interaction between the decorated particles 
on the conformation and distribution of the adsorbed polyelectrolyte chains, we evaluate the 
dipole moment of the complexes. It must be noted that for a system whose net charge is 
different from zero, the dipole moment is not univocally defined but assumes different values 
depending on the choice of the origin of the coordinates. However, to evaluate a 'degree of 
asymmetry' of the charge distribution of a pd-particle we calculate the expression 

f^='^qi{ri-rcc) (11) 

i 

where A''^ is the total number of charges g^ on the pd-coUoid and Tj individuate the position 
of each charge with respect to a fixed origin (in our case a corner of the box), and 

V^« I I ^ 
f _ 2^i W\^i (12) 

is the baricenter of the charges. 

Defining d = r^(.+ ~^cc- as the distance between the center of the positive charges r^j^ and 
that of negative charges r^^, (which, in our case, coincides with the center of the spherical 
macroion) Eq. [11] can be rewritten, after some algebra, as 

.--|f^<r (13) 

For an isolated pd-particle, where the charge distribution is, on the average, spherically 
symmetric, the dipole moment is zero since r^c+ = ^cc-- However, if the distribution of the 
adsorbed charges is 'distorted' by an external field, the modulus of fl increases, its value 
depending on the charge ratio and on the non zero value of d (see Eq. [T3l) . 

Due to the symmetry of the system, only the component of the dipole vector along the 
direction connecting the centers of the two macroions (the x axis with our choice of the 
coordinates) is affected. For 'intermediate' values of H {k,~^ < H < R^) depends on the 
net charge of the pd-particles. In Figure [7] /i^,. (if), fiy{H) and ^z{H) of two interacting 
complexes (A and B in Figure [9]) are shown for ^g = 0.26, Re = 5 nm k = 1 nm~^. 
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FIG. 7: The three components of the electrical dipole moment of the complexes A and B shown 
in Figure [9] as a function of the distance H between the macroions' surfaces; ^s = 0.26 . 

The increase of ^x-, for high values of H , is due to the asymmetry of the surface charge 
distribution induced by the opposing particle. At short distances a 'rapid' decrease of 
/ia; is observed, deriving from the entropic exclusion of the PE-chains from the contact 
region, where the steric interaction between monomers and between the monomers and the 
macroions' surfaces dominates. 

The effect of the charge ratio ^s on the dipole moment is shown in Figure O 
An asymmetric distribution of the mobile charges (the adsorbed polyelectrolyte chains) 
on the macroions can be generated only for distances between the complexes so short that 
the strength of the interactions is sufficient to promote energy driven diffusion processes 
that cause a redistribution of the adsorbed chains. The x-components of the dipoles are 
specular and their signs exchange in correspondence of the isoelectric point. In this way, 
the eventual aggregation of the two complexes occurs with the anti-parallel orientation 
of the dipoles (+- or -+ depending on wether the aggregation occurs below (^^ < 1) or 
above {^s > 1) the neutralization condition). This anti-parallel orientation of the dipole 
doublet signals the absence of a long-range dipole-dipole attractive contribution to the mean 
force inter-particle potential. Such findings give further support to the hypothesis that 
the attractive force that drives the aggregation of like charged pd-coUoidal particles arises 
from the charge mismatch between polyelectrolyte and polyelectrolyte-free domains on the 
opposing macroions' surfaces. 
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FIG. 8: A) the profiles of ficc{H) for botli tlie complexes (A,B Fig. [9]): D) ^s = 0.26 ; A ) ^s = 0.38; 
V) C. = 0.46; 0) 6 = 0.57; .) 6 = 1-0; ■) 6 = 1-2; T) 6 = 1.33. 
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FIG. 9: A snapshot of the labeled complexes (A,B) for ^s = 1-0 and the dipole orientation below 
and above the isoelectric point. 
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VI. FINAL REMARKS 

We studied the effective interaction and the induced charge asymmetry of PE-coUoid 
complexes within the Debye-Hiickel approximation for the case where the monomer size 
is small compared with that of the adsorbing macroion and the surface charge density of 
)are colloids is characteristic of the extensively investigated spherical DOTAP liposomes 



2|, ll9|, |25|, l42|-|45l. |58| . The interaction between non- neutral PE-decorated particles is purely 
repulsive for large screening length nR^ <^ 1 and it is well-described by a Debye-Hiickel 
potential. The interaction is sensibly modified as the screening is increased until a short 
ranged attractive component appears. 

From the competition between this short range attraction and the long range electrostatic 
repulsion due to the residual net charge of the like-charged complexes a global minimum in 
the pair potential arises whose position depends on the thickness of the adsorbed PE-layers. 
By adding chains on the macroions' surface the short range electrostatic component of the 
potential switches from repulsive to attractive well-below the neutralization condition (^g ~ 
0.46). For increasing values of the charge ratio ^si the adsorbed chains concomitantly reduces 
their typical dimension and form thick and spatially non homogeneous layers. Interestingly, 
the dipole moments of two approaching complexes are always anti-parallel (except at the 
neutralization point where all the dipoles' components are zero) so that the condensation of 
these complexes can not be attributed to dipole-dipole attractive interactions. 
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